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The problem of the equivalence of the spherical and mean spherical models, which has been 
thoroughly studied and understood in equilibrium, is considered anew from the dynamical point of 
view during the time evolution following a quench from above to below the critical temperature. 
It is found that there exists a crossover time t* ~ y 2 / d such that for t < t* the two models are 
^v^j ' equivalent, while for t > t* macroscopic discrepancies arise. The relation between the off equilibrium 

response function and the structure of the equilibrium state, which usually holds for phase ordering 
systems, is found to hold for the spherical model but not for the mean spherical one. The latter 
model offers an explicit example of a system which is not stochastically stable. 
PACS: 05.70.Ln, 75.40.Gb, 05.40.-a 
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I. INTRODUCTION 

Let us consider a system with scalar continous order parameter </?(x) in a volume V and hamiltonian 

n[y(x)] = \J dx{(V<p) 2 + r<p 2 (x)} (1) 



where r > 0. The spherical model of Berlin and Kac [1] is obtained by considering the extensive random variable 

o 
o 



# = i^dzV 2 (x) (2) 



and computing equilibrium properties with the gaussian weight p e [<p] — -^-e~T U ^ under the microcanonical con- 
straint ^ = a, where a is a given number. The mean spherical model of Lewis and Wannier [2], instead, is obtained by 
1^ ' imposing the constraint in the mean, or canonically ($) = a which makes the model considerably easier to solve. The 
spherical model was originally introduced by Berlin and Kac as an exactly soluble model displaying critical phenom- 
i ■ ena. Subsequently, the enforcement of a spherical constraint (in either form microcanonical or canonical) on the free 
part of non linear problems has become an extremely useful and practical way to generate mean field approximations. 

However, despite the great popularity of the method, it is usually overlooked that the equilibrium properties of 
the two models, as it was recognized early on, do coincide above but not below the critical point. The origin of the 
discrepancy was clarified first by Lax [3] and further investigated by Yan and Wannier [4] . Lastly, Kac and Thompson 
[5] showed how to connect averages in the two models. As it is easy to understand, the microcanonical and canonical 
constraint are equivalent as long as the fluctuations (S^ 2 ) = (^P 2 ) — (\E') 2 are negligible. Indeed, above the critical 
temperature Tc one has the usual behavior for thermodynamic quantities (S^ 2 ) ~ l/V. Not so below Tc, where 
q \ fluctuations of '5 turn out to be finite and independent of the volume (S 1 ^ 2 } ~ a 2 . 

An important consequence of this, as we shall see in the following, is that the nature of the mixed state below 
Tc is quite different in the two models. It is then interesting to investigate whether the two models are equivalent 
or not when considering the time dependent properties in the relaxation process following the quench from above 
to below the critical point. This is a relevant question since, in practice, dynamics can be solved only in the mean 
spherical case and in the literature it is given for granted that the mean spherical form of dynamics applies to the the 
spherical case as well. On the basis of the previous considerations, clarification of this point amounts to analyse the 
time evolution of (S^ 2 ). As we shall see, this depends on the order of the limits t — > oo and V — » oo. Namely, it turns 
out that if V is kept finite during the relaxation there is a crossover from the preasymptotic behavior (S^ 2 ) ~ l/V 
to the asymptotic one (5^ 2 ) ~ a 2 with the crossover time t* ~ y 2 / d . Instead, if V — + oo from the outset, then (S 1 ^ 2 ) 
stays negligible for any finite time since t* diverges. In any case, in the scaling regime t < t* , the time dependent 
evolution is identical in the two models, as usually assumed. 

Although reassuring, this conclusion opens an unexpected and interesting problem when it comes to test the connec- 
tion between static and dynamic properties introduced by Cugliandolo and Kurchan [6] for mean field spin glass and 
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then established in general by Franz, Mezard, Parisi and Peliti [7] for slowly relaxing systems. Dynamic quantities are 
the autocorrelation function C(t, t w ) and the integrated linear response function \(ty t w ), while equilibrium properties 
are encoded into the overlap probability function P(q) [8]. The statement is that if x{t,tw) depends on time through 
the autocorrelation function x(C(t,t w )), then one has 

where P(q) is the overlap probability function in the equilibrium state obtained when the perturbation giving rise to 
\{t,t w ) is switched off. Therefore, the unperturbed overlap probability P(q) can be recovered from dynamics if P(q) 
coincides or it is simply related to P{q). This happens for stochastically stable systems [7]. Formulated originally in 
the context of glassy systems, this connection between statics and dynamics applies also to phase ordering processes 
in non disordered systems [9]. A detailed account of the latter case can be found in Ref. [10]. Now, the interesting 
point is that the spherical and mean spherical model do share the same relaxation properties, therefore the same x(C) 
and the same P(q), while the corresponding unperturbed overlap probabilities are profoundly different. In particular, 
P(q) can be recovered from P(q) in the spherical model, but not in the mean spherical one. 

The paper is organized as follows. In Section 2 and in Section 3 the relevant properties of the spherical and mean 
spherical models are presented. The relation between the two models is discussed in Section 4 and comments on the 
connection between statics and dynamics are made in the concluding Section 5. 

II. SPHERICAL MODEL 



Dynamical evolution is described by the Langevin equation 

dy(x,t) 5H[<p] 



dt 5ip(x,t) 



n(x,t) (4) 



where r](x, t) is a gaussian white noise. Taking for H[<p] the gaussian model (1) and Fourier transforming with respect 
to space we have 

^tA = -u W {Kt)+r 1 (k,t) (5) 

with LOk = k 2 + r and 

MM)> = o 

( V (k, t) V (k', f )) = 2TV5 s+ 2 fi 5{t - t') W 
where T is the temperature of the quench. Integrating (5) we obtain 

tp(k,t) = R g (k,t)<p(k,0) + [ dt'R g (k,t-t')r){k,t') (7) 

Jo 

where R g (k,t) = e~ Wfct . Considering an infinite temperature initial state with 

MM)) = o 

Mfc,0Mfc',0)} = AV% +fe -, {) 

and taking averages over the initial state and thermal noise, from (7) follows (<p(k,t)} — 0, ((f(k,t)(p(k' ,t)) — 
C(k,t)V5^ + p where the equal time structure factor is given by 

C(k, t) = R 2 g (k, t)A + — [1 - R 2 g (k, t)}. (9) 

With the initial condition (8) and the linear equation (5) the configuration [ip(k,t)] executes a zero average gaussian 
process whose probability distribution is given by 



2 



Pg{[f(k)],t) = -l-e^&W^-'^M-S) (10 ) 
where Z g (t) = J %(fc)]e~^ 12s'P(S)C- 1 (U,tM-k) From ^ 

one can compute the one time properties of the gaussian 
model, including those at equilibrium obtained letting t — > oo. 

Let us next define the joint probability of a configuration [<p(k)] and of a random variable \& by 

Pg (Mk)]^,t) = Pg Mk)],t)s I *- -L 2 | V (fe) | 2 1 . (ii) 



Clearly, /9 S ([<£>(&)], i) is recovered integrating over while the probability of * is given by 

p g (9,t) = J d[<p(k)]p g {[<p(k)],9,t). (12) 
Introducing the probability of [<p(k)] conditioned to a given value of 

we may also write 

p g Mk)],t) = | d$p g (9,t)p g (* | [*>(*)],*)■ (14) 

Notice that conditioning with respect to <J/ is tantamount to imposing the spherical constraint. Hence, the probabilty 
distribution for the Berlin-Kac spherical model [1] can be written as 

p a ([<p(k)],t;a)=p g (9 = a\[<p(k)],t). (15) 



III. MEAN SPHERICAL MODEL 

As stated in the Introduction, the mean spherical model is obtained by imposing the constraint in the mean 
ytj J2k(\ I 2 ) = a - This can be done by the Lagrange multiplier method or, what is the same, by modifying (5) 
letting the parameter r to be a function of time 

^M = -[k 2 + r (t)Mk,t)+r 1 (k,t) (16) 
where r(t) is to be determined self consistently through the constraint, which can be rewritten as 

±£C(M) = a. (17) 

k 

The structure of the solution of the equation of motion remains the same 

<p(k, t) = R ms {k, t, 0)ip(k, 0) + / R ms (k, t, t')ri(k, t') (18) 

where now 



o 



R ms (k,t,t') = ^e- k2 ^ (19) 



with Y(t) — e®^ and Q(t) — J Q dt'r(t'). The equal time structure factor is given by 
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C(k, t) = R 2 ms (k, t, 0)A + 2T f dt'R 2 ms (k, t, t') 

Jo 



(20) 



where, clearly, the initial value A must also be consistent with (17). This requires V^^rA = 

I {2Tt) d exp(— fc 2 /A 2 )A = a where, as we shall always do in the following, in transforming the sum over k into 

an integral we make explicit the existence of an high momentum cutoff A. Hence, eventually A = (47r) d / 2 A~ d a:. 
Finally, the probability distribution keeps the gaussian form (10) 



Pms 



-2VE E v(fc)C'- 1 (MM-£) 



(21) 



where now C(k,t) is given by (20). 

In order to have an explicit solution the function Y(t) must be determined. This is done in Appendix I where, for 
simplicity, the computation has been limited to the case 2 < d < 4. For large time one finds 



Y 2 (t) 



b [e*/ r + (t/ry 



for T>T C 
for T < T c 



where 



2A d-2 



(47r) d / 2 (d-2) 



(22) 



(23) 



and the expressions of B,t,lo and t* are listed in Appendix I. Here we point out that B and r are independent of the 
volume for T > Tc, while B vanishes and r diverges as V —* oo for T < Tc- Inserting (22) in (20), for T > Tc we 
find 



C(k,t) 



^ e -2(fc 2 + *)t 

B 



T 



k 2 + h 



1 



_ 2(fe 2 + _L )t 



(24) 



and for T <T C 



C(k,t) 



A (t/t*) w e 



B [1 + [t/ryet/T] 



+ 2T / dt'e- 2k2 (t'/ty 



1 + {t'/t*Ye t '' T 
1 + (i/t*)"e*/ T 



(25) 



where t is the microscopic time necessary to elapse for (22) to apply. Notice that from Eq. (67) and from Eq. (74) 
follows that t* and r are both 0(V 2 ^ d ) for T = Tc, while for T < Tc and d > 2 the two time scales are separated 
with t* - V 2 ' d < t ~ V. 

In any case, r is the equilibration time. Taking t » r from Eq.s (24) and (25) follows 



C{k,t) = C eq {k) - 



fc 2 + r 



(26) 



where the equilibrium correlation length £ is related to t by £ 2 = 2r. Hence, using Eq. (67) Tc may be identified 
with the static transition temperature separating the high temperature phase where £ is independent of the volume, 
from the low temperature phase where £ diverges with the volume 




-l/(d-2) 



for < 2^ < 1 

for T = T c 
for T < T c . 



(27) 



Finally, let us comment on the nature of the equilibrium state. As Eq. (21) shows, each k mode is gaussianly 
distributed, with zero average, at any time including at equilibrium. Hence, in the low temperature phase the system 
does not order. The static transition consists in the k = mode developping a macroscopic variance as the temperature 
is lowered from above to below Tc [11] 



V for T > T c 
~ { for T = T C 

V 2 for T < Te- 



rn 
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For the sake of illustration let us consider T = where from C eq (k) = aVS^ Q follows 

1 2 

Pm.Mk)] = —==e-^ H 6(<p(k)). (29) 



V2naV 2 

fc#0 



IV. THE CONNECTION BETWEEN THE TWO MODELS 

In order to explore how the spherical and the mean spherical model are connected, let us rewrite Eq. (21) following 
the same steps which have led to Eq. (14) and obtaining 

p ms (&(k)},t) = J rf* Pms (*,t) Pms (* | [ip(k)],t) (30) 

which yields the relation between the two models since the conditional probability p ms (^ \ [ip(k)],t) enforces the 
constraint a la Berlin and Kac. Actually, to be precise, this quantity is not exactly the same as the spherical model 
distribution, since in Eq. (15) the constraint is imposed on the gaussian model, while here it is imposed on the mean 
spherical model. In the following we will ignore the difference. 

Looking at Eq. (30) the state in the mean spherical model can be regarded as the mixture of states in the spherical 
model with constraint values weighted by /9 ms ( x I , ,i). The properties of the two models are the same if this weight is 
narrowly peaked about the mean value (IS?) = a, while discrepancies are to be expected if the weight spreads over 
significantly different values of ^. Therefore, the key quantity controlling the connection between the two models is 
Pms^i). The corresponding characteristic function is 

G{x) = (e M *) = e-jEsMi-a-teC&t)] ^ 

where (•) stands for the average computed with (21). Moments of W then are given by (vp™) = d d /fJfJ ) and in 
particular 

<*(*)> = ^£C-(M) = « (32) 
fe 

(S^(t)) = (* 2 (t)) - (*(t)f = i^C 2 (M). (33) 

% 

Let us first consider what happens at equilibrium by letting t — > oo and inverting Eq. (31). Following Kac and 
Thompson [5] and using Eq.s (26) and (27) for T > T c one finds 

Pms (*) = -L= e -(*-«) 2 /2* (34) 

V 27T<7 

with a = J2k Ceq{k)- For large volume this can be rewritten as a = y J ^Jjd (jjrfgr=rp and since the integral is 
finite we have (S^ 2 ) ~ 0(1/V). Instead, for T <T C one finds 

for W < aT/Tc 

(tt-o,T/T c ) ,„_s 
e 2a(l-T/T c ) . [OO) 




for # > aT/Tc 



^2tx(^ -uT IT c )a(l-T IT C ) 

which gives 

(5V 2 ) =2[a(l-T/T c )} 2 . (36) 

Hence, we have that in going from the high temperature to the low temperature phase the fluctuations of \f from 
microscopic become macrocospic and, as anticipated above, significant differences must be expected in the equilibrium 
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states of the two models. In order to see this explicitely, let us consider T = where the equilibrium state in the 
mean spherical model is given by Eq. (29). The joint probability of [<p(k)] and is given by 

p ms ([<p(k)}, *) = -^= e-^ J] S (^))5 - 1 J dV(*)) (37) 



which, due to the presence of the 8 functions for k ^ 0, can be rewritten as 

fe#0 

Identifying the first factor in the right hand side with the T — > limit of Eq. (35), we find the conditional probability 



p ms (* | Mk)}) = 1 [% - ^V*) + <%o + 1] (39) 



k=£0 

which gives the T = state in the spherical model with a value \& of the constraint. Comparing Eq.s (29) and (39), the 
difference in the low temperature equilibrium states of the two models is quite evident. While in the mean spherical 
model ipo is gaussianly distributed, in the spherical model the probability of ip is bimodal. The large fluctuations (36) 
around the average of If do appear below Tq since it is necessary to spread out the weight of in order to reconstruct 
a gaussian distribution by mixing bimodal distributions. 

This difference in the structure of the low temperature equilibrium states shows up quite clearly in the corresponding 
overlap probability functions [8] 



P(q) = J d[<p]d[<P']p[<P]p[</]5(Q[<p, <p'] - q) (40) 

where 

Q[<P, f'] = ^ J d£<p(£)<p'W = ^ <p{k)<p?(-k). (41) 

k 

Keeping on considering, for simplicity, the T = states and looking at the characteristic function 6(A) = J dqP{q)e lXq 
from Eq.s (29) and (39) with * = a follows (Appendix II) 

<MA) = 1 (42) 
y/1 + (Xa) 2 

and 

9 S (A) = cos(Aa). (43) 

Inverting we have 

PmM = —K (\ q I /a) (44) 
an 

where K is the Bessel function of imaginary argument and 

P s (q) = l[8(q-a) + 8(q + a)}. (45) 

The plots of these two functions in Fig.l illustrate the great difference in the ground states. 

Having analysed the properties in the final equilibrium states, let us now go back to the dynamical problem. Notice 
that at t = 0, with C(k, 0) = A, from Eq. (31) follows that p ms (*,i = 0) is given by Eq. (34) with a = { l^ d % v ■ 
Hence, in a quench to T > Tq fluctuations of \& remain microscopic throughout the time evolution, while in a quench 
to T < Tc at some point (8^ 2 (t)) must cross over from 0(1/V) to 0(1). Again, for simplicity, we show this in the 
case of the T = quench. From Eq. (25) in the limit T^Owe have 
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C ^)= [i + wm (46) 



where to = (2A 2 ) -1 . Inserting into Eq. (33) we find for t/to > 1 and t* /to > 1 



which gives 



(S^(t))=2aHt/tr [ 0§0^ (47) 



(dV (t) > ~ ^ 2a 2 [l - (2 - 2- d / 2 )(t* /t) d l 2 ] for t » i*. (48j 



Therefore, t* is the crossover time separating the time regime t < f with (5^ 2 (t)) <~ 0(1/V) from the time regime 

* > t* with (a* 2 (t)) - o((r)- d / 2 ) - o(i). 

We may now comment on the non commutativity of the limits. If V is kept finite and the limit t — > oo is taken first, 
the system equilibrates and the two models at equilibrium are equivalent above Tq but not below. This conclusion 
keeps on remaining valid if, after having reached equilibrium, the limit V — > oo is taken. If, instead, the limit V — > oo 
is taken first then (d^> 2 ) ~ 1/V both for a quench above and below Tc, since in the latter case the crossover time t* 
diverges. Hence, the two infinite volume models are always equivalent during relaxation. Therefore, as stated in the 
Introduction, all dynamical quantities are the same in the two models if the V — > oo limit is taken from the outset. 



V. CONCLUDING REMARKS 



We have analysed the relationship between the spherical and the mean spherical model, at equilibrium and during 
relaxation from an initial high temperature state to a lower temperature state. By monitoring the behavior of the 
fluctuations of <]/ we have found that the two models are equivalent in a quench above Tc, while discrepancies arise in 
a quench below Tc if the volume is kept finite and the time is larger than the crossover time t* ~ V 2 / d . In the infinite 
volume limit t* diverges and the relaxation dynamics is the same in the two models for all times. In particular, the 
integrated auto-response function 



X(t,t w ) = 




is expected to have the same form in both models. As anticipated in the Introduction, this poses an interesting 
problem when considering the connection between statics and dynamics. 

In order to explain this it is necessary to expand somewhat on the behavior of the autocorrelation and response 
function in a phase ordering process [9,10,13]. The basic feature is the split of both these quantities into the sum of 
a stationary and an aging contribution 

C(t,t w ) = C st (t-t w ) + C ag (t/t w ) (50) 



x(t,t w )=Xst(t-t w )+x a g(t,t w ) (51) 

respectively due to thermal fluctuations and defect dynamics. The time scales of these contributions are widely 
separated. In particular C s t(t — t w ) decays rapidly from Mq — M 2 to zero (where M is the spontaneous magnetization 
at the temperature of the quench and Mo is the zero temperature spontaneous magnetization) while C ag (t/t w ) decays 
from M 2 to zero on a much longer timescale. The stationary terms in Eq.s (50) and (51) are related by the equilibrium 
fluctuation dissipation theorem 

T X st(t -t w ) = C st (t, t) - C st (t -t w ). (52) 

Rewriting the right hand side in terms of the full autocorrelation function one finds [9,10] 

f M 2 C(t, t w ) for M 2 <C< M 2 
1Xst[t ~ tw) -\M 2 -M 2 forC<M 2 [b6 > 

which yields 
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-^%P) c _ 9 = %-M 2 ). (54) 

The next statement is that the aging contribution of the response function obeys the scaling form 

Xag(t,t w )=t- a x(t/t w ) (55) 

and therefore vanishes for t w — > oo if a > 0. As of yet knowledge of the exponent a remains limited. According to 
euristic arguments [9,12] a ought to coincide with the exponent 9 controlling the defect density p(t) ~ t~ e , namely 
9 = 1/2 for scalar order parameter and 9 = 1 for vector order parameter [13], independently from dimensionality. 
However, exact analytical results for the one dimensional Ising model [14,15], careful numerical computations [10] for 
the Ising model in dimensions d = 2,3,4 and exact analytical results for the large N model [16] (which is equivalent 
to the mean spherical model) give the non trivial behavior 

f 0(-r=fy-) for d<du , ce x 
a = I y d v -d L > u 56 ) 

[6 for d > du y ' 

with logarithmic corrections at d = djj. Here d^ = 1 and djj = 3 for the Ising model, while dj, = 2 and djj = 4 in the 
large N case. What happens at d = dh, where a = 0, has been analysed in detail in [15,10,16]. Let us here consider 
d> dL 1 where a > and the aging contribution of the response function does asymptotically vanish. Then, putting 
together Eq.s (3) and (54) one finds [16,17] 

P(q)=5(q-M 2 ). (57) 

Now, as stated in the Introduction, if stochastic stability holds P(q) equals the overlap probability function P(q) of 
the unperturbed system, up to the effects of global symmetries which might be removed by the perturbation. This is 
what usually happens in a phase ordering system, like the Ising model, where the perturbation brakes the up down 
symmetry and one has 

P{q) = 26(q)P(q). (58) 

Even if only half of the states are kept in the order parameter function P(q) obtained from dynamics, using the 
symmetry P(q) = P(—q) it is obviously possible to recover the full unperturbed overlap function. Indeed, this is what 
takes place in the spherical model, where the unperturbed overlap function (45) is given by the sum of two 6 functions 
and P(q) = 29(q)P s (q) holds. Not so in the mean spherical case, where the unperturbed overlap function (44) is non 
trivial and P(q) ^ 29(q)P ms (q) as it can be seen at glance from Fig.l. Clearly, in the latter case P ms (q) cannot be 
reconstructed from knowledge of P(q). Therefore, stochastic stability does not hold in the mean spherical model. It 
might be interesting to investigate this point in other models treated with the spherical constraint. 

Acknowledgments - This work has been partially supported from the European TMR Network-Fractals c.n. FM- 
RXCT980183 and from MURST through PRIN-2000. 



VI. APPENDIX I 



In order to determine Y(t) explicitely, let us rewrite Eq. (17) by separating out of the sum the k = term, as we 
expect it to become macroscopically occupied at low temperature 



V 



For large volume the sum may be approximated by an integral and using Eq.s (19) and (20) we find 



aY 2 (t) - — dt'Y\t') 2T j o dt'f(t -t' + ^)Y 2 (t') = A 



where 



i , r d d k 

f{t+ 2^^IW) 



8*(*+2A*) 



-d/2 



(59) 



(60) 



(61) 
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The above equation can be solved by Laplace transformation obtaining 

A[l/V + zh(z)} 

C{Z) ~ az-2T[l/V + zh{z)] (62) 

where C(z) and h(z) are the Laplace transforms of Y 2 (t) and f(t + 5^7)1 respectively. The large time behavior of 
Y 2 (t) is controlled by the small z behavior of C(z). For 2 < d < 4 we have 

h(z) = K + 1 z d ' 2 - 1 + 0(z) (63) 
where K — (47r)~ d / 2 ^^- and 7 = — (87r)~ d / 2 r(l - d/2) are positive constants. Inserting in Eq. (62) we have 

A[l/V + Kz + ^z 4 ' 2 ] 



C(z) = 



(64) 



a(l - T/T c )z - 2T- f z d l 2 - 2T/V 
where Tc = ^ = (47r) d / 2 2^-2^ ■ Inverting the Laplace transform we have 

F 2 (t) = Be*/ r + /(t) (65) 
where the first contribution comes from the residue at the single pole at x = 1/r on the positive real axis and 

1 



I{t) = 



2iri 



dx[C(xe- l7T ) -£.(xe m )]e 



-xt 



(66) 



is the contribution from the cut along the negative real axis. Looking for the zero of the denominator of Eq. (62) we 
find 



2^(1 -T/T C ) 
( 7 y) 2 /d 

£{1-T/T C )V 



for < ^f- < 1 

for T = T c 
for T <T C 



(67) 



and computing the residue 



B = Aa {2Tr[a(d/2 - 1)(T/T C - 1) + dTr/V}} 1 



(68) 

For T > Tc the exponential dominates and the contribution from the cut can be neglected in Eq. (65), since to leading 
order B and r are independent of the volume. Not so for T < Tc, where taking into account the contribution from 
the cut we have 



with 



Y 2 (t) = B [e l ' T + {t/t*)- u [l ± (*/*)*]} 



( 2 -d/2 for T = T C 
\ d/2 for T <T C 



d/2 for T = T C 
1 for T < T c 



B = 



2 T-i 7 4/d for T ~T C 

^(l-T/Tc)- 2 for T<T C 



J_ (d/2-l)r 2 (2-d/2)V 
8n [2cos(d7r/2)r(2-d/2) 

a(i-r/r c )r(rf/2) ~ 

ATT(d/2-l)V 



2/d 



for T = T C 
for T <T C 



.* _ J \(d - 2) 2 A- 2 (87r) d - 2 7 8 / d ] 1/(d_4) V 2 ' d for T = T c 

for T < T c 



/ 2/d 
~87T 



(69) 

(70) 
(71) 

(72) 
(73) 
(74) 



and where in Eq. (69) the + and — sign apply to T = Tc and T < Tc, respectively. Notice that in all cases t* < r 
and t> t*, therefore the dominant contribution is given by 

t~ u for t < t* 
e*/ T for t>t*. 



Y 2 {t) 



(75) 
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VII. APPENDIX II 



Inserting Eq. (29) into the definition (40) we obtain for the overlap function in the mean spherical model at T = 

PmM = 2^ / dVod^e-^S^ - q). (76) 
The corresponding characteristic function is given by 

If 1 

©ms(A) = 2na y 2 / d<p (kp' cx P{"^72"(^o + <Pa - 2ia\ipov' )} (77) 

and going over to polar integration variables one finds 

6 ms (A) = 2iTa y 2 J J drr cxp{- (1 - 2iaX sin cos $)} 
1 f 277 d?9 

(78) 



7T 7 2 - aAe 4,? + aAe""? ' 

This integral can be rewritten in the form J- <£ - 22__ — _ where 7 is the circle of radius one with center in the 



origin of the complex plane. Since there is a simple pole at z = a\" inside 7 we obtain 

y/1 + (Xay 



and inverting the Fourier transform we find 

•'"-ss/-."^^- (80) 



ms\ 



Taking into account that there are branch points at ±i and that the integration contour is closed in the negative 
imaginary half plane for q > and viceversa for q < 0, eventually we obtain 

^ms(g) = — / dy^== = —K m) (81) 
where K is the Bessel function of imaginary argument. 
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